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Abstract 

At 6th order in perturbation theory, the random magnetic impurity problem at second 
order in impurity density narrows down to the evaluation of a single Feynman diagram 
with maximal impurity line crossing. This diagram can be rewritten as a sum of ordinary 
integrals and nested double integrals of products of the modified Bessel functions K^, and 

with v = 0,1. That sum, in turn, is shown to be a linear combination with rational 
coefHcients of (2^ — 1)C(5), J^uKQ{u)^du and u^Ko{u)^ du. Unlike what happens at 
lower orders, these two integrals are not hnear combinations with rational coefHcients of 
Euler sums, even though they appear in combination with C{5). On the other hand, any 
integral vT"^^ Kq{uY {uKi{u)Y du with weight p + q = & and an even n is shown to be a 
linear combination with rational coefficients of the above two integrals and 1, a result that 
can be easily generalized to any weight p + q — k. A matrix recurrence relation in n is built 
for such integrals. The initial conditions are such that the asymptotic behavior is determined 
by the smallest eigenvalue of the transition matrix. 



1 Introduction 



In Ref. [T], the quantum problem of a charged particle in a plane, coupled to a random Poisso- 
nian distribution of infinitely thin impenetrable Aharonov-Bohm flux tubes (magnetic vortices) 
perpendicular to the plane was considered. The relevant parameters are a = 4'/4>q, where (f) is 
the flux of a tube and 4)^ the flux quantum, and the mean impurity density p = N/V, where A'' 
is the number of impurities and V the area (in the thermodynamic limit, A^, ^ ^ oo with N/V 
fixed). Periodicity a ^ a + 1 and symmetry with respect to a = 1/2 imply that the A^-impurity 
partition function Z^r is invariant under a — 1 — a, and depends only on a(l — a). 
One was interested in the average partition function 
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where Zq is the free partition function. 

With account for the a dependence of Zjy, Eq. ^ becomes a sum of terms proportional 
to p"'Q^ with m > n. For small a, the leading terms are (/oa)"- They yield the partition 
function of the charge in the mean magnetic field pacpQ, which replaces the local magnetic field 
(pYliLi ^{^ — ^i), being the location of the i-th impurity. Terms with m > n are perturbative 
corrections to the mean-field expansion, which originate from disorder effects. For the 1-impurity 
case p], which is exactly solvable, one obtains (Zi) — Zq = a{a — l)/2. For the 2-impurity case, 
nontrivial Feynman diagrams with maximal impurity line crossing appear at order p^oi^, i.e., an 
electron interacting 4 times with 2 impurities; at order /J^a^, i.e., an electron interacting 6 times 
with 2 impurities, etc. Knowing the p^a^ Feynman diagram is sufficient to get a rather precise 
analytical estimate of the critical disorder coupling constant Uc — 0.35, above which oscillations 
in the density of states, corresponding to the Landau levels in the mean magnetic field picture, 
disappear. That is a clear indication that the system becomes fully disordered. Note, on the 
other hand, that at weak disorder, when a < Qc, the broadening of the Landau levels due to 
disorder fits nicely into the Integer Quantum Hall Effect paradigm. 

In Ref. [I], the p^a^ diagram was reduced to a multiple temperature integral 

V„.= TdA t\lh f\pri-.r,^^^^m^) (3) 
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where a = Pi — (32, b = (32 — (3^, c = (3^ — (3/^, d = (5^^ — (3i + (3. A direct step-by-step integration 
gave 

,^^.^0>(L.m\ (4) 
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that is, a linear combination of 1 and C(3) = 7C(3)/2 with rational coefficients. One inferred by 
connecting the integral ^ to the partition function that 
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Thus, an Euler sum of level 3 emerges, which fits into the general scheme of Feynman diagram 
expansion in perturbative field theory [3], where Euler sums are known to play a central role. 
These sums are defined as 

C{ni,n2,...,np) = 

where n = ni + n2 + ■ ■ ■ + Up is the level of the sum. At level n, the simplest sums are 
C(n) = EZi and Ca{n) = EZi ("1)'/^", with (2" - l)C(n)/2- = (C(n) - Ca{n))/2. 
On the other hand, Eqs. ©-((II) imply that the nontrivial part of /p2Q,4 is 
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In Ref. [4j, algebraic manipulations and a Laplace transform with respect to /5 let one factorize 
this multiple integral as 

-/ dadbdcdd dt — e-('^+*+'=+'^)-*(o + 6 + c + d) = / uKQiuf{uKi{u))'^du 

2 Ja,b,c,d=0 Jo Ca Jq 
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where Ky{u) are the modified Bessel functions {t = /2). Consequently, 
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More generally, it was shown in by means of a simple integration by parts that any integral 
of the form 

[•oo 

u''+^KQ{uy{uKi{u)ydu (7) 





with weighJll p + q = A and n even, is a linear combination with rational coefficients of C(3) and 
1. This was achieved via a 2 x 2 matrix recurrence relation for the integrals u^~^^KQ{u)'^(iu 
and u^^^ Ki{u)'^du, with n > 4 (see Sec. 3 for a derivation of this recurrence and its straight- 
forward generalization to any weight k = p + q). A remarkable result is the fact that the initial 
condition 

5 4 -27 + 7C(3) /""5^.^4, 53 - 9C(3) 
u^'Koiu) du = — , / u^'Kiiu) du= — (8) 

Jo 
happens to be such that the asymptotic behavior of the recurrence relation for large n is governed 
by the smaller of the two eigenvalues {1/16, 1/4} of the asymptotic recurrence matrix. 

In Ref. [5], these considerations were extended to integrals involving Ki, as well as Ii, — 
specifically, the set ^ where either one of the Kq^s is replaced by Iq or one of the i^i's is 
replaced by Ji (so one has a product of three K^'s and one 1^). These integrals were shown, 
again via an elementary integration by parts, to be linear combinations with rational coefficients 
of 3C(2) and 1. The 2x2 recurrence matrix is identical to the one obtained in the previous case, 
up to a minus sign in the off-diagonal elements. The initial condition 
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u^Ko{ufloiu)du=^, ru^K,iufh{u)du = ^ (9) 



is such that the asymptotic behavior is governed, as expected, by the bigger of the eigenvalues 
{1/16, 1/4} — it does not for sure coincide with the unique initial condition jS]) associated with 
the smallest eigenvalue. 



2 The p^a^ diagram with maximal impurity hne crossing 

Integrals of products of modified Bessel functions appear to play a central role in the perturbative 
analysis of the 2-impurity problem. Indeed, whereas the p^a^ diagrams can be easily shown to 
reduce to p'^a^ diagrams, the relevant p'^a^ diagram with maximal impurity line crossing is much 
more arduous to compute. Following the same route as for the p^a^ diagram, and again taking 
a Laplace transform, one has obtained \E] the expression 
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■^Here and in the sequel we define the weight of an integral of a product of Bessel functions K,^ as the total 
power of the K^^s. 
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It contains not only, in analogy with order p^a^, the integral 



uKQ{u)^{uKi{u)fdu 

of weight 6, but also a particular combination of nested double integrals of products of modified 
Bessel functions and ly. Components of those double integrals, if integrated indivually from 
to oo, are, as already demonstratecil, linear combinations of either ^(3) and 1 or of 3^(2) and 
1. For example, in addition to Eq. jH]), one has 

{uKi{u))h{u)KQ{ufdu = ^ 

o 

uKQ{u){uK^{u)fh{u)du = ^^^^ 
°^ u'KQ{u)\uK^{u))h{u)du = ^^1^ 

-J^ 

In effect, there is a mapping, via an integral, of a product of — or a product of Ky and 
—onto a linear combination with rational coefficients of C(3) and 1 — or C,{2) and 1, respectively: 
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For the double integrals, the same scheme is at work, but now the mapping is to a "polyzeta" 
object: 

/•oo i^u 

if,9)^ / fiu)du / gix)dx = C[f,9] 



Since 

/•oo 

/(n)dn / g{x)dx = j f{u)du I g{x)dx — I g{u)du I f{x)dx 



one has 

C[/,5] = C[/]Cb]-Cb,/] 

in analogy with the relation involving the standard polyzeta function ({p, q) = X]n>m nPm^' 

C(p, q) = C{p)Ciq) - dp + q)- aq,p) (11) 

By analogy with lower orders, one might expect the p^a^ corrections to be a linear com- 
bination with rational coefficients of Euler sums up to a certain level. The structure of the 
double integral in Eq. (fTOl) clearly indicates that the highest level should be 5. Indeed, the 
constituent single integrals in the first and second terms reduce to levels 3 and 2, respectively; 
by virtue of Eq. (fTT]l . the product C{3)({2) is associated with C(5)- The last term in ([TO]) . 

uKo{u)^{uKi{u))^ du, is a Bessel integral of weight 6. Given that an integral ([6|) of weight 
4 is related to C(3), i-e., level 3, this agains hints at level 5 in the case at hand. 

However, a search for an integer relation with the PSLQ algorithm does not confirm this 
expectation. On the contrary, it indicates that 

1 1 31 

W = ^y^ ^^o(n)^dn+-y^ u^Ko{ufdu-—Ci5) (12) 



^Note, however, that x Ko{x)'^Ki{x)^ dx would be divergent if one set it = 0. 
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is a linear combination with rational coefficients of not only, as expected, a level 5 Euler sum 
(2^ — but also of two numbers of weight 6, uKo{u)^ du and u^Kq{u)^ du, neither 

of which is a linear combination with rational coefficients of Euler sums of level 5. 

Moreover, the same PSLQ search shows that the last term in Eq. (fTOll is a linear combination 
of the same two numbers: 



2 f'°° 1 
uKo{u)^{uKi{u)fdu = — / uKoiufdu-- 
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(13) 



Hence, also the sum of double integrals in (iTOl) is by itself a linear combination with rational 
coefficients of the same three numbers that appear in (fT2]l — specifically, — uKq{u)^ dn/10 + 
Ko{uf du/4 - 31C(5)/160. 



3 Integrals u''+^Ko{u)P{uKi{u)y du with even n 

Equation (fT3l) for weight 6, together with the integrals of weight 4 in [HIS], are but particular 
cases of a much more general pattern involving integrals of the form u'^'^^ Ko{uy {uKi{u)y du 
with even n. 

At a given weight k = p + q, denote 
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where j = 0,1,2, k. For the integral to be finite, n > j — 1 is required. Integration by 
parts, taking into account that dKQ{u)/du = —Ki{u) and d{uKi{u)) / du = —uKq{u), gives a 
recurrence relation 
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One has dei An^ 



■ /^i n\} \"-!'\} \ "'t^'' 1 I ^^ when k is odd, and when k is even, meaning that 

(n+2){n+l)(n)...(n—k+2) ' ' o 



in the latter case the -^^^ s are linearly related. Indeed, for an even k, one has 



k/2 . , 

j:y)\n-2i+2)(^fyi%=o 



(17) 



We choose here, as a matter of convenience, to start the recurrence at n = fc. 
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By applying relation (fTSj) twice, one obtains a mapping {-f„^2 ji ~^ ^-^n j}' 

(fc) ^ (n + l)(n + 2)(j-l)j (,.) 

(n-j + 2)(n- j+4) "+2.i-2 



(n + l)(n + 2) 
n - j + 2 



r(fc) 



n — j + 2 n — j + 4 

, (n + l)(n + 2)(fc-j-l)(fc-i) (,) 

+ (n-j + 2)2 ^"+2j+2 U«j 

This can be inverted, as long as takes into account the linear relation (fT71) in the case of an even 
k. 

The recurrence relations (fTSl) . (fTSl) conserve the parity of n — j; thus, all In j^s are divided 
into two subsets, and the relations operate separately within each subset. For reasons that will 
become clear later, we focus on one of those — the one for which n — j is even. To span this 
subset, it is enough to assume that n is even, then take the set 

Infl^ ^2' ^nl (^™) (19) 

[remembering that these integrals are linearly related via Eq. (fTTl) ]. or 

i'l I^l^ n>k-l (fcodd) (20) 

and note that the integrals 

^+i,3> (fceven) (21) 

and 

^+i,3> (^odd) (22) 

are related to ([H]) and ([20]) . respectively, by inverting the recurrence ([T5]l [again, taking into 
account ([T71) if /c even]. The union of sets (fT9]l - (l22]) . at a given k, is tantamount to the family of 
integrals 

u''+^KQ{uY{uKi{u)fdu (23) 

with weight p + q = k and an even n. 

For an even k, by virtue of Eq. ([TS]) . all integrals from the family ([T9]l . ([2T]) can be expressed 
as a linear combination with rational coefficients of the initial conditions 

fik) Ak) j{k) ■^ 

-'fc.O' -'A:,2' ^k,k l^^^ 

still taking into account Eq. ((TT]). Noting further that Eq. ([15]) imphes 21^^^-^^ -,^ = 

(A; — l)-^^*2_2) and that, trivially, = 1/^', one can trade /^'^ for 1. Using ([T7]) . one can 

get rid of one more element: thus, all integrals ([T9]) . ([2T]) are linear combinations with rational 
coefficients of a basis made of the k/2 independent (in the sense that none of them is a linear 
combination with rational coefficients of the others) numbers 

n jik) jik) Ak) . .2rN 
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Likewise, for k odd, starting with the initial conditions 

Ak) Ak) Ak) i^r,-. 

^k-1,0^ ^fc-1,2' • • • ' ^k-l,k-l y-^^l 

and again using /^^^^ = one concludes that any integral in ([20]) . ([22]) is a linear combi- 

nation with rational coefficients of the basis made of the {k + l)/2 independent numbers 

n Ak) Ak) Ak) . ^27^ 

Last but not least, one can also show, by applying (fTSj) appropriately for < n < /c, that for 
k even, the basis (^5]) can be mapped on the basis 

ri Ak) Ak) j{k) . 

t-*^ ' -'0,0 ' -'2,0 ' • • • ' -'fe-4,0-1 '.^"z 

and for k odd, the basis ([27ll can mapped on 

r-i 7-(fe) Ak) Ak) . 

i-*- ' -'0,0 ' -'2,0 ' • • • ' -'fc-3,oJ y'^^) 

Therefore, any integral in the set (f23ll is a linear combination with rational coefficients of the 
basis (f28l) or ([29]) . for A; even or odd, respectively. 

Consider now the asymptotic regime, n oo. In that limit, Eq. (fTSl) becomes 



n,j — ^j,j-2^n+2,j-2 ^n+2,j ^i,j+2-'n+2,j+2 V'J'-'j 

where the elements of the tridiagonal transition matrix A^^^ (indexed so that their subscripts are 
always even) no longer depend on n: 

^55-2 = (j - l)i > ^ = 2j{k-j) + k, Af]^^ = {k-j-l){k-j) (31) 

The [{k — k mod 2)/2 + 1] eigenvalues of this matrix are: A;^, {k — 2)^, {k — 4)^, . . . ; the last one 
is 1 for odd k and for even k. In the latter case, before inverting the matrix, one has to reduce 
its dimension by one, using Eq. (|17p . Thereupon, the largest eigenvalue of the inverse matrix is 
1 or 1/4 for /c odd or even, respectively, whereas the smallest one is l/k"^. Remarkably, as we 
discovered experimentally, with the initial conditions (f24l) . (f26]l it is the smallest eigenvalue that 

(k) 

determines the asymptotic behavior of /^^ j. In the asymptotic regime itself, this can be under- 
stood by noting that the eigenvector corresponding to said eigenvalue, i.e., to the /c^ eigenvalue 
of A^^\ is {1, 1, . . . , 1} — because Af^._^ + Af^- + Af^-^^ = k^. Now, in the limit n ^ oo, the 

integral I^^j does not depend on j, because the integrand u"~^^ Ko{u)''~^ Ki{uy peaks at large 
values of u, where both Kq{u) and Ki{u) can be approximated by their common asymptotic 
behavior, K^{u) — ^ \/^^~"- Therefore, the vector {1^ ]} becomes, in the asymptotic limit. 



2u^ • -^ii^i^i^i^; "ii^ v^^L^i n,j ^ 

proportional to {1, !,...,!}. The fact that the initial condition ([2^ . ([26]) leads to this special 
asymptotic behavior, as well as the fact that said condition can be expressed in terms of the 
basis means that the building blocks of (pSIl - iP^ . namely Bessel function integrals 

u"'~^^ Kq{u)^ du with n even, may play some special role in number theory — like Euler sums 
do (but again, these integrals are not rational linear combinations of those sums). 

As an example, we detail the recurrence relations for weight /c = 3, where one has computed 

m 



(3) 



0,0 



1 1 2 



7 



One has, for an even n > 2, 



n,0 



(3) 



/ 3(ra + 1 
n + 2 



Vn,2/ 



6(n + 1) \ 



n + 2 

2(n + l) (n + l)(7n + 6) 



The inverse relation is 
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which in the asymptotic Hmit turns into 
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9 3/ 
with eigenvalues {1/9, 1}. 

Clearly, as alluded to above, instead of the set (f23ll . one could have focused on the set 
u"^^Ko('u)^(tii^i('u))'' du with weight p + q = k and n odd — corresponding to the subset of 
integrals (fT4l) with odd n — j. The recurrence still operates within this set, and the same kind 
of algebra as above is at work. However, the integrals in this set 

• do not play any role in the perturbative analysis of the random magnetic impurity problem 
(at least up to 6th order); 

• but still lead to an asymptotic behavior governed by the smallest eigenvalue of the cor- 
responding asymptotic recurrence matrix. This can be easily seen in the case A; = 3, 
where the initial conditions /g^j, I'^fl for the integrals J^^g and J^^j with n even and the 
corresponding recurrence relation lead to an asymptotic governed by the smallest of the 
eigenvalues {1/9,1} of the asymptotic matrix. 

(n+l) 



Note finally that at weight A; = n + 1, by singling out the simplest integral /q q 
(l28l) or (l29]l. one arrives at the number 



1 



K[n] 



(n + l)! 



in the basis 



(36) 



with k(0) = 1; k(1) = 1/4; k{2) = 45/3!, see Eq. (I32|); k(3) = 7C(3)/(8 x 4!), etc. As already 
said, this number is analogous to Euler sums of level n, but it is not a rational linear combination 
of those. When n ^ 00, the l/(n + 1)! normalization in (f36l) is such thailll 



g-2r'(i) Yi^ai n{n) = lim (C(n) - 1) 



(37) 



When n —1, on the other hand, n{n) = 1/(1 + n)^ plus logarithmic subleading terms and a 
constant. Clearly, the function K,{n) defined for n > — 1 real can be analytically continued to the 
function k{s) defined on the whole complex half-plane K(s) > — 1. 



'^This can be easily shown by recognizing that when n - 
origin, where Kq(u) can be approximated by Ko{u) ~ — log 



00, the integrand uKo{u)"''^^ has a peak near the 
It follows that the main contribution to 



the integral is u iiro(w)"^^ dit ~ Jg 
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""•"Mu, which trivially yields JST 
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4 Conclusion 



We have demonstrated that a class of integrals involving Bessel functions, which arise in perturba- 
tion theory — in particular, in the two-dimensional problem of random magnetic impurities — can 
be expressed, via recurrence relations, as linear combinations with rational coefficients of "basis" 
integrals. Some of the latter, in turn, reduce to Euler sums, but most do not. Additionally, 
these same basis integrals, for an even power of the argument, turn out to generate the unique 
initial conditions for the recurrence relations which make the smallest eigenvalue of the transi- 
tion matrix determine the asymptotic behavior. The nature of this phenomenon has yet to be 
understood more deeply. 

It is not only single but also some double nested integrals that happen to be linear combina- 
tions with rational coefficients of said basis integrals: one that does is the integral that figures 
in the 6th order of perturbation theory in the physical problem at hand. Understanding what 
other double integrals fall into the same class remains another open question. 

Acknowledgements: One of us (S.O.) would like to thank Jean Desbois for discussions and 
some technical help, in particlular at the end of Section 3. 

Note added: After completion of this work, we became aware of Ref. [8] and in particular 
of Ref. 26 therein, where some results overlap with those of section 3. 
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